Strong dineutron correlation in 8 He and 18 C 
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We study the spatial structure of four valence neutrons in the ground state of 8 He and 18 C nuclei 
using a core+4n model. For this purpose, we employ a density-dependent contact interaction among 
the valence neutrons, and solve the five-body Hamiltonian in the Hartree-Fock-Bogoliubov (HFB) 
approximation. We show that two neutrons with the coupled spin of S=0 exhibit a strong dineutron 
correlation around the surface of these nuclei, whereas the correlation between the two dineutrons 
is much weaker. Our calculation indicates that the probability of the (IP3/2) 4 and [(IP3/2) 2 (pi/2) 2 ] 
configurations in the ground state wave function of 8 He nucleus is 34.9% and 23.7%, respectively. 
This is consistent with the recent experimental finding with the 8 He(p, t) 6 He reaction, that is, the 
ground state wave function of 8 He deviates significantly from the pure (IP3/2) 4 structure. 
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I. INTRODUCTION 

It has been well recognized that the pairing correla- 
tion and couplings to the continuum spectra play an es- 
sential role in weakly bound nuclei [1, 2]. Although the 
dineutron structure as a consequence of the pairing cor- 
relation has been suggested for some time in 11 Li and 
6 He nuclei[l, 3], it is only recently that a strong indi- 
cation of its existence has been obtained experimentally 
in the Coulomb dissociation of 11 Li [4]. The new mea- 
surement has stimulated lots of theoretical discussions on 
the dineutron correlation, not only in the 2n halo nuclei, 
11 Li and 6 He [5-8], but also in medium-heavy neutron- 
rich nuclei [9, 10] as well as in infinite neutron matter 
[11, 12]. 

In Ref. [7], we have studied the behaviour of valence 
neutrons in 11 Li at various positions from the center to 
the surface of the nucleus. We have found that i) the two- 
neutron wave function oscillates near the center whereas 
it becomes similar to that for a bound state around the 
nuclear surface, and ii) the local pair coherence length has 
a well pronounced minimum around the nuclear surface. 
This result clearly indicates that a strong di-neutron cor- 
relation between the valence neutrons is present on the 
surface of the nucleus. 

An important next question is how the spatial struc- 
ture of valence neutrons evolves from that in the 2n-halo 
nucleus, 11 Li, when there are more numbers of neutrons. 
Although Refs. [9, 10] have partially addressed this ques- 
tion by studying a two-particle density for medium-heavy 
neutron-rich nuclei, one would also need to explore a 
four-particle density, or many-particle density in general, 
in order to shed light on the ground state properties of 
neutron-skin nuclei. 

For this purpose, 8 He makes the most suitable nucleus 
to study. 8 He is expected to have the a+4n structure 
[13-20], and thus provides a bridge between the 2n-halo 
nuclei and heavier skin nuclei. We mention that the spa- 
tial structure of the four valence neutrons in 8 He has been 
discussed in Ref. [16], where the authors constructed the 



ground state wave function by assuming that the four 
neutrons occupy the IP3/2 state in a harmonic oscillator 
potential. However, this model is too simplistic, since it 
completely neglects the pairing correlation and the con- 
tinuum couplings. Notice that the mixing of many partial 
wave components, especially those with different parities, 
is essential in order to have a spatially compact dineutron 
structure [9, 10, 21]. In fact, we do not see any indica- 
tion of dineutron correlation in the result of Ref. [16] (see 
Fig. 2 of Ref. [16]), despite that the dineutron structure 
is expected to be enhanced in many neutron-rich nuclei 
[7, 9, 10]. 



The purpose of this paper is to reinvestigate the spatial 
structure of the four valence neutrons in 8 He by taking 
into account consistently the pairing and the continuum 
effects. To this end, we use the core+4?i model, and di- 
agonalize the five-body Hamiltonian in the Hartree-Fock- 
Bogoliubov (HFB) approximation. We also study the 18 C 
nucleus as another nucleus which is expected to have the 
core+4n structure[22-25]. We will demonstrate below 
that the pairing correlation leads to the strong dineutron 
structure in 8 He and 18 C, in contrast to the result of Ref. 
[16]. 



The paper is organized as follows. In Sec. II, we de- 
tail the HFB method based on the core+4n model. In 
Sec. Ill, we apply the method to 6 He, where the result 
of the exact diagonalization of the three-body (a+n+n) 
Hamiltonian has been obtained [6, 26]. We compare the 
HFB result with the exact result, and discuss the appli- 
cability of the HFB method for the study of the spatial 
structure of valence neutrons. In Sec. IV, we present the 
results for the 8 He and 18 C nuclei. We discuss the two- 
and four-particle densities, as well as the probability of 
the single-particle components in the ground state wave 
function. We then summarize the paper in Sec. V. 



2 



II. HFB METHOD FOR A CORE+4™ MODEL 

A. HFB equations 

In order to study the structure of the 8 He and 18 C 
nuclei, we employ the corc+4n model, and consider the 
following Hamiltonian: 
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Here, m^v is the nucleon mass, A is the mass number of 
the nucleus, V n c is a potential between a valence neutron 
and the core nucleus, and v nn is the pairing interaction 
among the valence neutrons. T cm is the kinetic energy 
for the center of mass motion of the whole nucleus. In 
Refs. [6, 26], the center of mass motion is treated exactly 
by introducing the recoil kinetic energy of the core nu- 
cleus (see Ref. [14] for the derivation of the recoil term). 
In this paper, we approximate the treatment by taking 
only the diagonal components in T cm , as is often done in 
mean-field calculations [27] (notice that the off-diagonal 
components contribute only to the exchange part of the 
mean- field potential). This leads to the renormalization 
of the nucleon mass, m = Aj \A — 1) • ttin- 

Although the five-body Hamiltonian (3) could be di- 
agonalized exactly e.g., with the stochastic variational 
method [15], we seek an approximate solution using the 
HFB method [2, 28-30]. The ground state wave function 
in the HFB method is given by [28] 



|HFB)=IJ&|0), 



where the quasi-particle operator (3k is given by 



fa = 
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In this equation, er = ±1/2 is the spin coordinate, a\. a is 
the creation operator of nucleons at the position r and a, 
and Uk and Vk are the HFB quasi-particle wave functions. 
In this paper, we employ a density-dependent pairing in- 
teraction [1] for v nn given by 



,(r,r') = V„ 1 
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Po 
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where r = (r + r')/2 and p*(r) = pc( r ) + Pv( r ) is the 
total density, pc and p v being the density of the core 



nucleus and the valence neutrons, respectively. For this 
interaction, the expectation value of the Hamiltonian (3) 
with the HFB state (4) reads [2] 



E = (HFB | .ff | HFB), 



J dr (J^ T ( r ) +V n c(r)p v (r 
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where the kinetic energy density r(r), the particle density 
p v (r), and the pairing density p v (r) are given by 

< r ) = EEi w *( r ' ff )i 2 ' ( 9 ) 
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respectively. In deriving Eq. (8) , we have used the prop- 
erties of time-reversal symmetry [2] . 

The equations for the HFB wave functions Vk and Uk 
are obtained by taking the variation of the energy ex- 
pectation value (8) with respect to the particle and the 
pairing densities. This leads to the HFB equations in the 
coordinate space representation [2, 29, 30], 
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where A is the Fermi energy. The mean-field Hamiltonian 
h is given by 
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while the pairing potential A(r) is given by 
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We solve the HFB equations (12) sclf-consistently by 
expanding the HFB wave functions on the eigen functions 
of the mean-field Hamiltonian h [31]. In doing so, we 
respect the Pauli principle and explicitly exclude those 
states which are occupied by the neutrons in the core 
nucleus. Notice that the HFB method could be applied 
to light neutron-rich nuclei without introducing the core 
nucleus [32, 33]. We nevertheless treat only the valence 
neutrons explicitly, since it is not straightforward to sep- 
arate between the core and the valence parts from the 
HFB ground state wave function (4). 
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B. Two- and four-particle densities 

In order to discuss the spatial structure of the valence 
neutrons, we compute the two- and four-particle densities 
using the solution of the HFB equations. Using Wick's 
theorem, the two-particle density can be expressed as [9] 

p 2 (ra,r'a') = <HFB|4 ff 4 v , arV 'ar CT |HFB}, (17) 
= \p v (ra,r'a')\ 2 -\p v (ra,r'a')\ 2 
+p v (ra)p v {r'a'), (18) 

where a = —a, and the off-diagonal components of the 
densities are given as (see Eqs. (10) and (11)), 

p v (ra,r'a') = £ V k (r, a)V k *(r', a'), (19) 
k 

Pv (ra,r'a') = -£V fc (r,(7)I/ fc *(rV). (20) 
k 

In order to evaluate the four-particle density, 

P4(xi,X2,X 3 ,Xi) 

= (HFB | a Xi a\ 2 a\ 3 a Xi a Xi a X3 a X2 a Xl | HFB) , (21) 

where x = (r,a), we find it useful to express the HFB 
ground state wave function, Eq. (4), using the canonical 
basis. The canonical basis function ip is the eigenfunction 
of the density matrix (19) and satisfies [2, 28] 

I dr' Pv (ra, rV) ^ P (r', a*) = v 2 p %{r, a). (22) 

a' J 

In this paper, we construct the canonical basis by ex- 
panding ip p on the HF basis, as is done for the HFB 
wave functions (see the previous subsection). Using the 
canonical basis, the HFB ground state wave function is 
given in the BCS form as [2, 28], 



with 
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where u p = ^Jl — v 2 p and p is the time-reversed state of 

p. Since the creation operator for the canonical basis, 
a p , is related to the creation operator in the coordinate 



space, a\, a , as 



a l= dr Yl ^p( r > CT ) a TV> ( 25 ) 
Eq. (24) is transformed to [2], 

HFB) cxcxp J drdr' ^2z(r<r,r' a') a\. a al„ a ^j 



10), 
(26) 



Z(ra, r'a>) = -2 £ ^ ^(r, a)^ p (r', a 1 ). (27) 
p>0 p 



To evaluate the four-particle density, (21), we first per- 
form the particle number projection onto the HFB state, 

fV|HFB) cx {-^ J drdr'J2 z (rT,r'a')al. a al., a ^j |0). 

(28) 

The four-particle density is then obtained as 

p 4 (x 1 ,x 2 ,x 3 ,x 4 ) cx \f(x 1 ,x 2 ,x 3 ,x 4 )\ 2 , (29) 

with 

f(x 1 ,x 2 ,x 3 ,x 4 ) 
= (Z(x 1 ,x 2 ) - Z{x2,x\)){Z{x 3 ,Xi) - Z(x 4 ,x 3 )) 
+ (Z(xi,x 3 ) - Z(x 3 ,Xi))(Z(x4,x 2 ) - Z(x 2 ,x 4 )) 
+ (Z( Xl ,x 4 ) - Z(x i ,x 1 ))(Z{x 2 ,x 3 ) - Z(x 3 ,x 2 )\30) 



C. Probability for shell model configurations 

Using the canonical basis representation of the HFB 
state, (23), one can also calculate the probability for 
a shell model configuration, [(kk)(k'k')], for the four- 
particle systems when k and k' represent the canonical 
basis states. It is given by 

P kVc „ = -l|((fcfc)(fc'fc')|HFB)| 2 = ±v 2 v 2 , J] u 2 p , 

P ^k.k'(>0) 

(31) 

where the normalization factor M reads 
M = (HFB 1^=4 1 HFB) = -L (" d4e-«+ T\(u 2 p +e 2 ^v 2 

(32) 

Here, we have used the explicit form of the number pro- 
jection operator [28], 
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If the angular momentum components are explicitly 
expressed, the probability (31) reads 

1 /ft 



P ^=Jf\ 2 )(*&)W'- a II (^'^',(34) 



for the (lj) 4 configuration, while 



x J! («?»i")°'"> ( 35 ) 

l"j"jilj,l'j> 

for the (lj) 2 (l'j') 2 configuration with lj ^ l'j'. In these 
equations, ftj = (2j + l)/2 is the pair degeneracy for the 
angular momentum j state. 
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III. APPLICABILITY OF HFB METHOD: 
ON 6 HE NUCLEUS 



TEST 



Let us now numerically solve the HFB equations and 
discuss the spatial structure of neutron-rich nuclei. Be- 
fore we do this, however, we first examine the applicabil- 
ity of the HFB method by applying it to the three-body 
model of 6 He nucleus [6, 26]. This model has been solved 
exactly by diagonalizing the Hamiltonian matrix. A com- 
parison of the HFB solution with the exact result for this 
model will provide an idea on whether the HFB method 
is good enough to discuss the dineutron correlation in 
neutron-rich nuclei. 

To this end, we use the same neutron-core potential, 
V n c, and the same pairing interaction v nn as in Refs. 
[6, 26]. The neutron-core potential is taken as a Woods- 
Saxon form. Since we use the renormalized mass m (see 
Eq. (3)) instead of the reduced mass, we multiply the 
factor (A — 1)/A ■ Ac /{Ac + 1), where Ac is the mass 
number of the core nucleus, following the prescription 
given in Ref. [26] . The pairing interaction is given as 



v nn (r, r') = 5(r - r') (v a + - 



+ cxp[(r - R p )/a p 

(36) 

We multiply an overall scaling factor to this interaction 
so that the two-neutron separation energy of 6 He is repro- 
duced with the HFB method. We use the same value for 
all the other parameters as in Ref. [6] . Note that the last 
term in Eq. (14) disappears for the pairing interaction 
given by Eq. (36), since the interaction does not depend 
explicitly on the density, but the density dependence is 
parameterized by the Fermi function. 

Figure 1 shows the two-particle density for 6 He in the 
5 = channel, that is, pii^i ],r 2 !). As we have done 
in Ref. [6], we set n = r 2 = r and plot the density as 
a function of r and the relative angle between the spin 
up and down neutrons, 9. Figure 2 shows the same two- 
particle density, but we multiply the factor 87r 2 r 4 sin# 
[6]. The top panels in these figures show the exact so- 
lution of the three-body Hamiltonian [6] , while the mid- 
dle panels are for the HFB results. One can clearly see 
that the localization of the two-particle density around 
9 ~ in the three-body model is well reproduced by the 
HFB method, although the HFB density has a somewhat 
longer tail and the density around 9 ~ n is largely sup- 
pressed. The localization of the two-particle density is 
nothing but the manifestation of the strong di-neutron 
correlation in a halo nucleus 6 He. The longer tail of the 
HFB density may be due to the asymptotic behavior of 
the pair density p v (ri 1,r 2 j) in Eq. (18), which is dif- 
ferent from that of the normal density p v {r\ t) T "2 I) [2]. 
The similarity between the exact and the HFB results 
for the two-particle density is rather striking, and it is 
clear that the HFB method can be utilized to discuss, 
at least qualitatively, the strong dineutron correlation in 
neutron-rich nuclei. 

We also study the effect of particle number projection 
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FIG. 1: (Color online) The S = component of the two- 
particle density for 6 He as a function of ri = T2 = r and the 
angle between the valence neutrons, 6. The top panel shows 
the exact solution of the three-body model, while the middle 
panel is obtained with the HFB method. The bottom panel 
shows the result of the HFB + particle number projection. 



on the two-particle density. The bottom panels in Figs. 
1 and 2 are obtained by taking the number projection 
onto the HFB ground state (with the variation before 
projection (VBP) scheme [34]) in a similar way as in Eq. 
(28). The two-particle density thus obtained is not nor- 
malized and the scale is different between the middle and 
the bottom panels. However, we can see that the depen- 
dence of the two-particle density on r and 9 is almost the 
same between the two panels. Therefore, we conclude 
that the effect of number projection on the two-particle 
density is rather small as far as the two-particle density is 
concerned, although the projection might still affect the 
density if the variation after projection (VAP) scheme is 
employed. 

Table 1 summarizes the occupation probabilities for 
the 6 He nucleus. Although the absolute value is some- 
what smaller, the HFB well reproduces the dominance of 
the (P3/2) 2 configuration in the ground state wave func- 
tion. Again, the HFB method provides a good estimate 
of the ground state properties of neutron-rich nuclei even 
when the number of particle is as small as two. 
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FIG. 2: (Color online) Same as Fig.l, but with a multiplica- 
tive factor of 87r 2 r 4 sin0. 

TABLE I: Comparison of the exact and the HFB results for 
the occupation probabilities in the ground state of 6 He. 



configuration 


exact [6] 


HFB 


(si/% 


2 


3.04 % 


7.25 % 


(Pl/2 


) 2 


4.85 % 


9.21 % 


(P3/3 


) 2 


83.0 % 


57.4 % 


(ds/a 


) 2 


1.47 % 


3.86 % 


( d 5/2 


) 2 


6.11 % 


6.85 % 


(f5/2) 


2 


0.035 % 


2.30 % 


(f 7 /a) 


2 


0.075 % 


3.33 % 



IV. DINEUTRON CORRELATION IN HE AND 

18 C 

We now solve the HFB equations for the 8 He and 18 C 
nuclei. We use the same neutron-core potential, V n c, for 
8 He as in Refs. [6, 26] , while we use the set D in Refs. [24, 
25] for the 18 C nucleus. As in the previous section, these 
potentials are scaled by a factor of (A — l)/A- Ac / (Ac + 
1). For the core density, pc, we use those in Refs. [16, 
35] . We determine the strength of the pairing interaction 



so that the experimental ground state energy relative to 
the corc+4n threshold, E = -3.112 MeV for 8 He and 
-10.385 MeV for 18 C, is reproduced with p =0.32 fm" 3 
(i.e., the mixed pairing interaction[36, 37]). With the 
cut-off energy of e cu t + A=40 MeV in the single-particle 
space, this leads to Vq = —502 MeV fm 3 for 8 He and 
V = -538 MeV fm 3 for 18 C. 



TABLE II: The results of the HFB calculation for the Fermi 
energy A and the root-mean-square (rms) radius, r rms , for the 
s He and 18 C nuclei. 



nucleus 


E g . a . (MeV) A (MeV) 


r rms (fm) 


r& (fm) 


8 He 


-3.112 -0.0715 


3.23 


2.49 ± 0.04 [13] 


is c 


-10.514 -2.522 


2.92 


2.90 ± 0.19 [38] 



The results of the HFB calculation are summarized in 
Table II. Although our purpose in this paper is not to 
reproduce the experimental data, but to discuss qualita- 
tively the dineutron correlation in 8 He and 18 C, we notice 
that the root-mean-square radius for the 18 C nucleus is 
well reproduced with the present calculation. The mean- 
field potential in Eq. (14) is shown in Fig. 3, in which 
the dashed and the solid lines correspond to the neutron- 
core potential, V n c, and the total mean-field potential, 
respectively. The difference between the two potentials 
originates from the effect of pairing correlations among 
the valence particles on the mean field potential. As a 
consequence, the mean-field potential for 8 He posses one 
bound single-particle state while the neutron- 4 He poten- 
tial V„c alone does not hold any bound single-particle 
state, reflecting the Borromean nature of the 6 He nu- 
cleus [6, 26]. For the 18 C nucleus, the same effect shifts 
the single-particle energy from —1.072 to —1.768 MeV 
for the 2si/ 2 state and from -0.414 to -1.664 MeV for 
ld 5 / 2 - 

TABLE III: Probability of a few shell model configurations in 
the HFB ground state wave function for 8 He and 18 C. 



nucleus 


configuration 


probability (%) 


8 He 


[(1P3/2) 4 ] 


34.9 




[(l P 3/2) 2 (Pl/2) 2 ] 


23.7 




[(P3/2) 2 (d 5 /2) 2 ] 


10.7 




[(s 1 /2) 2 (P3/2) 2 ] 


7.8 


18 c 


[(ld 5 /2) 4 ] 


32.2 




[(ld 5/2 ) 2 (2s 1/2 ) 2 ] 


26.2 




[(ld 5 /2) 2 (d 3 /2) 2 ] 


11.8 




[(d 5/ 2) 2 (f7/ 2 ) 2 ] 


7.17 



The probability for a few single-particle components 
in the ground state wave function is listed in Table III. 
For the 8 He nucleus, although the largest probability is 
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FIG. 3: The mean-field potential for 8 He (the upper panel) 
and for 18 C (the lower panel). The dashed line shows the 
neutron-core potential, V„c, while the solid line is for the total 
mean-field potential in the solution of the HFB equations. 



found in the [(IP3/2) 4 ] configuration, the other compo- 
nents also have an appreciable probability. Therefore, 
this nucleus largely deviates from the pure [(IP3/2) 4 ] con- 
figuration, in accordance with the recent experimental 
finding with the 8 He(p,i) 6 He reaction [39]. For the 18 C 
nucleus, the ground state wave function mainly consists 
of the [(ld 5 / 2 ) 4 ] and the [(ld 5/ / 2 ) 2 (2si/ 2 ) 2 ] configurations, 
while the [(ld 5/2 ) 2 (d 3/2 ) 2 ] and the [(d 5/2 ) 2 (f 7/2 ) 2 ] config- 
urations are also appreciable in Table III. 

The top panel in Figs. 4 and 5 shows the two-particle 
density, p 2 (r, r — 0, f;r, r, j), for 8 He and 18 C, respec- 
tively. The middle panels in these figures show the same 
two-particle density, but with the multiplicative factor 
of 87r 2 r 4 sin#. For both the 8 He and 18 C nuclei, one 
clearly finds a strong concentration of two-particle den- 
sity around 9 ~ at around the nuclear surface. This is 
similar to what has been found in the Borromean nuclei, 
11 Li and 6 He [6, 26] (see also Figs. 1 and 2), and indi- 
cates clearly the strong dineutron correlation [7, 9, 10] in 
these nuclei. 

Since the strong dineutron structure is apparent for a 
spin-up and a spin-down neutrons in these nuclei, we next 
plot the four-particle density for the two-dineutron con- 
figuration, that is, the four-particle density with x\ = 
(r,f = 0, f), x 2 = (r,f = 0,|), x 3 = (r,r,|), and 
X4 = (r,r, l) in Eq. (21). This is plotted in the bot- 
tom panels in Figs. 4 and 5 for 8 He and 18 C, respec- 
tively. For the 8 He nucleus, the four-particle density for 
the dineutron-dineutron configuration has a peak around 
9 ~ 7r/2. A similar result has been obtained with a three- 
body model calculation with the dineutron clusters, that 
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FIG. 4: (Color online) The two-particle density, p2(r,r — 0, T 
; r, r, J.), for the 8 He nucleus as a function of ri = T2 = r and 
the relative angle 6 between a spin-up and a spin-down neu- 
trons (the top panel). The middle panel shows the same two- 
particle density multiplied by a factor 8n 2 r 4 sin 9, while the 
bottom panel is for the four-particle density for the dineutron- 
dineutron configuration. 



is, a + n 2 + n 2 [17]. The peak around 9 ~ 7r/2 arises 
from the main component of the wave function, that is, 
the [(lp 3 / 2 ) 4 ] configuration, for which the four-particle 
density is proportional to sin 4 9 oc | Yn | 4 . For the 18 C nu- 
cleus, the four-particle density has two peaks, one around 
9 ~ 54 deg. and the other around 9 ~ 118 degree. This 
can again be understood in terms of the [(ld 5 / 2 ) 4 ] con- 
figuration, for which the four-particle density is propor- 
tional to (3|r 22 | 2 + 2|r 21 | 2 ) 2 . 

To demonstrate more clearly the similarity between the 
four-particle density to that for the main components, the 
bottom panels of Figs. 6 and 7 show the four-particle 
density for the [(lp3/ 2 ) 4 ] and [(ld 5 / 2 ) 4 ] configurations in 
the neutron-core potential V n c, respectively. To this end, 
we adjust the depth of the neutron-core potential so that 
the energy of the lp 3 / 2 and ld 5 / 2 states is a quarter the 
energy of 8 He and 18 C, respectively. The similarity be- 
tween the four-particle density for the correlated wave 
functions (Figs. 4 and 5) and that for the uncorrelated 
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FIG. 5: (Color online) Same as Fig. 4, but for 18 C. 



wave functions (Figs. 6 and 7) is apparent. This is a 
natural consequence of a short range nature of nuclear 
interaction: the two neutrons with the same spin have 
to be far apart in space due to the Pauli principle and 
thus their distance is likely larger than the range of the 
nuclear interaction. As a consequence, the interaction 
between the two dineutrons becomes weak, despite that 
the two neutrons in the same dineutron having different 
spins strongly interact with each other. From this consid- 
eration, we conclude that the two dineutrons are moving 
rather freely in the core+4n nuclei respecting solely the 
Pauli principle. 

Note that the pairing interaction yet plays an essen- 
tial role in the two-particle density. Without the pairing 
correlation, the two-particle density for the uncorrelated 
wave functions has symmetric bumps both around 9^0 
and 9 ~ 7r, as is shown in the upper panel in Figs. 6 and 
7 (see also the middle panels, that show the two-particle 
density with the weight of 87r 2 r 4 sin 9) . The pairing corre- 
lation mixes several angular momentum components in 
the ground state wave function, eliminating essentially 
the bump around 9 ~ it. The configuration mixing of 
different parity states seen in Table III is essential to 
have the di- neutron peak in the middle panel of Figs. 4 
and 5. 
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FIG. 6: (Color online) Same as Fig. 4, but for the uncorrelated 
[(lp 3/ / 2 ) 4 ] configuration for the 8 He nucleus. 



Another way to investigate the four-particle density is 
to plot the density distribution of the second dineutron 
when the first dineutron is put on the z-axis [9], rather 
than assuming that the distance from the core nucleus 
is the same between the two dineutrons. The top, mid- 
dle, and bottom panels of Fig. 8 show the four-particle 
density of 8 He for the dineutron-dineutron configuration 
when the first dineutron is at z—l.b, 2.5, and 3.5 fm, re- 
spectively The pairing correlation is taken into account 
in the plot. The same plot for the 18 C nucleus is shown 
in Fig. 9. These figures demonstrate that the distance of 
the second dineutron from the core, r2, increases as the 
distance of the first dineutron, n, increases, tending to 
7"i ~ r%. The angular distribution of the second dineu- 
tron, on the other hand, is almost independent of the 
position of the first dineutron. This behaviour is consis- 
tent with the four-particle density shown in the bottom 
panels of Figs. 4 and 5. 



V. SUMMARY 

We have discussed the dineutron structure in the 2n 
halo nucleus 6 He as well as in the core+4n nuclei, 8 He 
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FIG. 7: (Color online) Same as Fig. 5, but for the uncorrelated 
[(ld 5/ / 2 ) 4 ] configuration for the 1S C nucleus. 



and 18 C. For this purpose, we employed the density- 
dependent contact interaction among the valence neu- 
trons, and diagonalized the core+xn Hamiltonian (x—2 
for 6 He and x — 4 for 8 He and 18 C) with the Hartree- 
Fock-Bogoliubov (HFB) method. From the comparison 
with the exact solution of the three-body Hamiltonian for 
the 6 He nucleus, we found that the HFB method is sat- 
isfactory enough to discuss the spatial structure of the 
valence neutrons. For the 8 He and 18 C nuclei, we in- 
vestigated both the two- and the four-particle densities. 
We showed that the two-particle density takes the largest 
value when the spin-up and the spin-down neutrons are at 
the same position, that is nothing but the manifestation 
of the strong dineutron correlation. With this result in 
mind, we particularly discussed the four-particle density 
for the dineutron-dineutron configuration. We found that 
two dineutrons weakly interact with each other, simply 
respecting the Pauli principle. The four-particle density 
of the HFB calculation in fact resembles to that for the 
uncorrelated wave functions. This result is entirely due 
to a short range nature of nuclear interaction. Namely, 
the two neutrons with the same spin have to be far apart 
in space due to the Pauli principle and thus their distance 
is likely larger than the range of the nuclear interaction. 
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FIG. 8: (Color online) The four-particle density of 8 He for the 
dineutron-dineutron configuration when the first dineutron is 
on the 2-axis. The top, middle, and bottom panels correspond 
to the cases where the first dineutron is at z=1.5, 2.5, and 3.5 
fm, respectively. 



As a consequence, the interaction between the two dineu- 
trons becomes weak, while the two neutrons in the same 
dineutron strongly interact with each other. 

We have also discussed the probability for the single- 
particle configurations in the ground state wave func- 
tion. Our HFB calculations indicate that the 8 He nu- 
cleus consists of the [(IP3/2) 4 ] configuration by 34.9% 
and of the [(lp3/2) 2 (pi/2) 2 ] configuration by 23.7%, 
while the 18 C nucleus consists of the [(ld 5 / 2 ) 4 ] and the 
[(ld 5 /2) 2 (2si /2 ) 2 ] configurations by 32.2% and 26.2%, re- 
spectively. The result for the 8 He nucleus is consis- 
tent with the recent experimental finding with the two- 
neutron transfer reaction, 8 He(p, i) 6 He, that indicates 
an appreciable mixture of the configurations other than 
[(IP3/2) 4 ], e.g., [(lp 3 /2) 2 (Pi/2) 2 ], in the ground state of 
8 He. It would be interesting to analyse the experimental 
data for the 8 He(p, t) 6 He reaction with the wave function 
obtained in this paper. This will be a topic for a future 
publication. 
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FIG. 9: (Color online) Same as Fig. 8, but for 18 C. 



We thank J. Dobaczewski for useful discussions. This 
work was supported by the Japanese Ministry of Educa- 
tion, Culture, Sports, Science and Technology by Grant- 
in- Aid for Scientific Research under the program number 
19740115. 



[1] G.F. Bertsch and H. Esbensen, Ann. Phys. (N.Y.) 209, [10] 
327 (1991). 

[2] J. Dobaczewski, W. Nazarewicz, T.R. Werner, J.F. [11 
Berger, C.R. Chinn, and J. Decharge, Phys. Rev. C53, [12 
2809 (1996). 

[3] M.V. Zhukov, B.V. Danilin, D.V. Fedorov, J.M. Bang, [13 
I.J. Thompson, and J.S. Vaagen, Phys. Rep. 231, 151 
(1993). ' [14 

[4] T. Nakamura, A.M. Vinodkumar, T. Sugimoto, N. Aoi, 

H. Baba, D. Bazin, N. Fukuda, T. Gomi, H. Hasegawa, N. [15 
Imai, M. Ishihara, T. Kobayashi, Y. Kondo, T. Kubo, M. 
Miura, T. Motobayashi, H. Otsu, A. Saito, H. Sakurai, S. 
Shimoura, K. Watanabe, Y.X. Watanabe, T. Yakushiji, [16 
Y. Yanagisawa, and K. Yoneda, Phys. Rev. Lett. 96, 
252502 (2006). [17 

[5] F. Barranco, P.F. Bortignon, R.A. Broglia, G. Colo, and 

E. Vigezzi, Eur. Phys. J. All, 385 (2001). [18 

[6] K. Hagino and H. Sagawa, Phys. Rev. C72, 044321 

(2005); C75, 021301 (2007). [19 

[7] K. Hagino, H. Sagawa, J. Carbonell, and P. Schuck, Phys. [20 
Rev. Lett. 99, 022506 (2007). 

[8] C.A. Bertulani and M.S. Hussein, Phys. Rev. C76, [21 
051602 (2007). 

[9] M. Matsuo, K. Mizuyama and Y. Serizawa, Phys. Rev. [22 
C71, 064326 (2005). 



N. Pillet, N. Sandulescu, and P. Schuck, Phys. Rev. C76, 
024310 (2007). 

M. Matsuo, Phys. Rev. C73, 044309 (2006). 

J. Margueron, H. Sagawa, and K. Hagino, Phys. Rev. 

C76, 064316 (2007). 

I. Tanihata, D. Hirata, T. Kobayashi, S. Shimoura, K. 
Sugimoto and H. Toki, Phys. Lett. B289, 261 (1992). 
Y. Suzuki and K. Ikeda, Phys. Rev. C38, 410 (1988); Y. 
Suzuki and Wang Jing Ju, Phys. Rev. C41, 736 (1990). 
K. Varga, Y. Suzuki, and Y. Ohbayasi, Phys. Rev. C50, 
189 (1994); K. Varga, Y. Suzuki, and R.G. Lovas, Nucl. 
Phys. A571, 447 (1994). 

M.V. Zhukov, A.A. Korsheninnikov, and M.H. Smed- 
berg, Phys. Rev. C50, Rl (1994). 

A.V. Nesterov, V.S. Vasilevsky, and O.F. Chernov, Phys. 
Atom. Nucl. 64, 1409 (2001). arXiv:nucl-th/0006001. 
A. Adahchour and P. Descouvemont, Phys. Lett. B639, 
447 (2006). 

Y. Kanada-En'yo, Phys. Rev. C76, 736 (2007). 
N. Itagaki, M. Ito, K. Arai, and S. Aoyama, to be pub- 
lished. 

F. Catara, A. Insolia, E. Maglione, and A. Vitturi, Phys. 
Rev. C29, 1091 (1984). 

Y. Suzuki, H. Matsumura, and B. Abu-Ibrahim, Phys. 
Rev. C70, 051302(R) (2004). 



[23] W. Horiuchi and Y. Suzuki, Phys. Rev. C73, 037304 

(2006) ; ibid. 74, 019901 (2006). 

[24] K. Hagino and H. Sagawa, Phys. Rev. C75, 021301(R) 

(2007) . 

[25] N. Vinh Mau, Nucl. Phys. A592, 33 (1995). 

[26] H. Esbensen, G.F. Bertsch and K. Hencken, Phys. Rev. 

C56, 3054 (1999). 
[27] M. Bender, K. Ruth, P.-G. Reinhard, and J.A. Maruhn, 

Euro. Phys. J. A7, 467 (2000). 
[28] P. Ring and P. Schuck, The Nuclear Many Body Problem 

( Springer- Verlag, New York, 1980). 
[29] J. Dobaczewski, H. Flocard, and J. Treiner, Nucl. Phys. 

A422, 103 (1984). 
[30] A. Bulgac, e-print: nucl-th/9907088. 
[31] K. Hagino and H. Sagawa, Phys. Rev. C71, 044302 

(2005). 

[32] J. Meng and P. Ring, Phys. Rev. Lett. 77, 3963 (1996). 
[33] M.V. Stoitsov, J. Dobaczewski, W. Nazarewicz, S. Pit- 

tel, and D.J. Dean, Phys. Rev. C68, 054312 (2003); 

http:/ /www.fuw.edu.pl/ dobaczew/thodri/thodri.html. 



10 



[34] K. Hagino, P.-G. Reinhard, and G.F. Bertsch, Phys. Rev. 

C65, 064320 (2002). 
[35] Yu.L. Parfenova, M.V. Zhukov, and J.S. Vaagen, Phys. 

Rev. C62, 044602 (2000). 
[36] J. Dobaczewski and W. Nazarewicz, Prog. Theor. Phys. 

Suppl. 146, 70 (2002). 
[37] J. Dobaczewski, W. Nazarewicz, and M.V. Stoitsov, 

Euro. Phys. J. A15, 21 (2002). 
[38] E. Liatard, J.F. Bruandet, F. Glasser, S. Kox, T.U. Chan, 

G.J. Costa, C. Heitz, Y. Elmasri, F. Hanappe, R. Bim- 

bot, D. Guillemaudmueller, and A.C. Mueller, Europhys. 

Lett. 13, 401 (1990). 
[39] N. Keeley, F. Skaza, V. Lapoux, N. Alamanos, F. Auger, 

D. Beaumel, E. Becheva, Y. Blumenfeld, F. Delaunay, A. 

Drouart, A. Gillbert, L. Giot, K.W. Kemper, L. Nalpas, 

A. Pakou, E.C. Pollacco, R. Raabe, P. Roussel-Chomaz, 

K. Rusek, J.-A. Scarpaci, J.-L. Sida, S. Stepantsov, and 

R. Wolski, Phys. Lett. B646, 222 (2007). 



